Independent Noise Approximation for Spin-Boson Decoherence 
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Quantum error correction is a solution to preserve the fidelity of quantum information encoded 
in physical systems subject to noise. However, unfavorable correlated errors could be induced even 
for non- interacting qubits through the environment (bath), when they are "packed" together. The 
question is, to what extent can we treat the noise induced by the bath as independent? In the 
context of the spin-boson model, we show that, under some reasonable constraints, the independent 
noise approximation could be valid. On the other hand, in the strongly correlated limit, we show 
how the method of decoherence free subspace can be made applicable. Combining these two methods 
makes fault-tolerant quantum computation promising in fighting against correlated errors. 

PACS numbers: 03.67.Pp, 03.65.Yz, 03.67.Lx 



Quantum information encoded in realistic physical sys- 
tems is very fragile, compared with the classical counter- 
parts. Quantum error correction [1, 2] is possible and ef- 
ficient provided that the errors are independent or weakly 
correlated. For independent noise, where errors are mod- 
eled as stochastic events, error thresholds [3, 4, 5] can be 
estimated to give an upper bound of the error rate for 
fault-tolerant quantum computation. The effects of cor- 
related errors are generally not favorable for error cor- 
rection as more restrictive error thresholds are needed. 
In reality, quantum noise may not behave like stochastic 
noise — the independent noise assumption is question- 
able. No physical qubit can be genuinely isolated from 
the environment (bath). In particular, for solid state sys- 
tems, if the qubits are interacting with the same bath, 
indirect correlations could be generated, and the protec- 
tion from error correction codes becomes ineffective. Our 
goal here is to illustrate the essential features and gen- 
eral structures of the environment induced decoherence 
for the purpose of quantum error correction. 

This work is motivated by the recent studies [4, 5, 6, 
7, 8] related to this problem. However, some physical 
considerations have been ignored, and consequently the 
adverse effects of the bath on quantum error correction 
were often overestimated. In [4, 5, 6], attempts of deriv- 
ing mathematically rigorous error thresholds were made 
by summing over the so-called "fault paths". However, 
for general environments, the spin-bath coupling H-ffssH 
is unbounded and hence fails the analyses. The spin- 
boson model [9] (with [H s ,Hsb] = cf. Eq. (3)) have 
been considered in [7, 8], where a coarse-grained descrip- 
tion on the dephasing effects were studied in [7] , but the 
possibility of reducing the collective effects (also ignored 
in [8]) through decoherence free encoding were not con- 
sidered. Moreover, the effects of including non-trivial lo- 
cal Hamiltonians Hs have not been properly explored, as 
we shall see it does play an important role in generating 
spatial error correlation. This work is aimed at provid- 
ing a more complete and more physical picture on this 
problem, and suggests a more optimistic view on the ef- 



fects of environment induced decoherence. Particularly, 
we shall address two questions: (1) To what extent, can 
the noise induced by the bath be considered as indepen- 
dent? (2) For fully correlated noise, how do we minimize 
the decoherence effects? 

We formulate the problem in terms of the relevant two- 
point correlation functions. For a generic boson bath 
(e.g. acoustic phonon bath) and physical qubits with 
local energy splittings A, we found that the error cor- 
relation generated from the energy-converving processes 
(here called "bit-flipping" case) becomes significant when 
the qubit-qubit separations are less than a length scale 
A* = Tic/ A, where c is the wave speed of the bath. The 
physical origin of the spatial error correlation is due to 
the constructive interference of the local disturbances 
from the qubits to the bath. This provides an answer 
to question (1). For question (2), we found that when 
A — > (here called dephasing case), the errors are fully 
correlated when the long wavelength modes of the bath 
become dominant. Provided that the effective qubit- 
qubit interactions [cf. Eq.(7)] are properly handled, it 
is possible to minimize the effects by using the noiseless 
subspace (or decoherence free subspace DFS) to encode 
the qubits. 

Definition of Independent Noise — We first clarify the 
concept of error correlation in the context of quantum 
error correction. The degree of error correlation may be 
best quantified by the deviation from the independent 
noise. We motivate the definition of 'independent noise' 
through a simple example: consider a register of qubits, 
consisting of N spin- 1/2 particles, initialized in some en- 
coded state \Qubits), subject to a random classical mag- 
netic field which causes an unwanted evolution 
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where Wk = {Xk,Yk, Z^} is one of the Pauli matrices 
(Wf = I) acting on spin k and the random variable 
|efc| <§; 1 is assumed to be small. In powers of each 
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term in this series is in general not mutually orthogonal. 
However, after error syndrome detection, the final state 
will be projected into some states with definite parity 
with respect to the corresponding stabilizers. The prob- 
ability P\(Wk) of single error on spin k can be estimated 
by the amplitude square Ai(Wk) 2 = ||efcWfc \Qubits)\\ = 
|efc| 2 of the Wk term. The amplitude of any two-qubit 
error is simply a product of the amplitudes of two single- 
qubit errors A 2 [W 3 W k ) = \e 3 e k \ = A x {W j )A 1 (W k ), for 
anv j 7^ ^ Provided that all |efc| <C 1, the standard 
threshold theorem about independent stochastic noise 
should be valid. To generalize the above argument to 
the errors caused by the noise from the environment (for- 
mally denoted as \Env)), we may (naively) promote the 
complex numbers e k to operators e k acting on the envi- 
ronment. Then, Ai(W k ) 2 = (Env\ e\e k \Env) and so on. 
However, generally A 2 (WjW k ) ^ A 1 {W j )A 1 (W k ), but 
we want to know under what circumstances the equality 
would hold. This motivates our definition of the inde- 
pendent noise model: 

A n (WiW 2 ...W„) « A 1 (W 1 )A 1 (W 2 ) ■ ■ ■ A 1 (W n ) (2) 

for all possible combinations of the Pauli matrices. The 
sign offers us the flexibility to neglect higher order 
corrections, if any. Below we shall study the question 
about the validity of the independent noise approxima- 
tion in the context of the spin-boson model. The method 
to be introduced may also be applicable to more general 
models of the environment, and shall likely give the same 
qualitative results. 

Spin-Boson Model — We assume that the interaction 
between the boson bath and the physical qubits is suf- 
ficiently "weak" to be considered as perturbation. The 
full Hamiltonian H = Hs + Hb + Hsb, or 

consists of three parts, respectively, Hs (to be speci- 
fied later) the system Hamiltonian of N physical qubits, 
Hb the bosonic bath Hamiltonian, and Hsb the inter- 
action Hamiltonian, Z k = X^=i Zje ' rj , where fj is 
the position of the j spin. For simplicity we assume that 
all of the spins are coupled uniformly g 3 = g with the 
bath. Morevoer, we shall confine our attention to the 
baths having two generic properties: (a) the modes are 
3-dimcnsional k = (k x ,k y ,k z ) and (b) the modes have a 
linear dispersion relation u> k = ck. 

Case I: Dephasing — To continue with a more quan- 
titative analysis, here we shall first consider the Hs = 
case, since it is exactly solvable [12, 13] (alternatively, 
take the A — > limit from Eq. (14)) and contains 
rich physical contents. The unitary operator (K = 1) 
U = cxp(-iHt) can be decomposed into three parts 



where Hb is defined in Eq. (3), 
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G SB (t) = Y d Z j 4> j {t) , 

3=1 

and with f k ( rj ,t) ee (g/uj k )e-^ (l - e*"**), 
<t>j (*) = E [? k ( r J'*) a l ~ h ( r i>*)* ak 



(5) 



(6) 



The effective spin-spin interaction is of the ZZ form: 

,2 



H eff = E ^2^ k Zl [uj k t - sin (uj k t)] 



(7) 



It is generated through the virtual transitions of the spin- 
bath interaction. By definition, it does not depend on 
the state of the bath and the effects are predictable, and 
hence should not be considered as noise [10]. In principle, 
it can be included as part of the quantum circuit, or elim- 
inated through some refocusing schemes well-developed 
in NMR quantum computing. For the moment, we shall 
neglect it. The term cxjp(-iHBt), containing no spin 
variables, is irrelevant for evaluating the error amplitudes 
and will also be neglected. Our goal is to expand the re- 
maining parts in a power series, similar to Eq. (1), of 
g (or Zj), and then estimated the relative size of each 
term. 

We assume the modes are isotropic, i.e., cu- k = uj k . 
Then all commutators [</>j{t),<f)k(t)] = 0, containing odd 
parity terms sin [A: ■ (fj — r k )], should vanish. This means 



cxpGss(i) = '[[exp[Z j <t>j(t)] 



(8) 



which resembles Eq. (1). Suppose at t = the spins are 
in a pure state and the bath is in a thermal state. The 
square of the error amplitudes are given by (<f)' = —<f>) 



A n (ZiZ2---Z n )~ 



(9) 



From the Bloch identity [11], we have (exp£j<f)j) = 
exp [(£j£ m /2) (4>j4> m )]- Here repeated indices imply sum- 
mation, and £j is an arbitrary complex number. The 
quantity A\ can then be evaluated by the standard tech- 
nique of summing different "pairing" of the two-point 
correction functions, which are the key quantities to be 
considered for error correlation, 

((f> 3 (f> m ) =]T [h k (R,t)(ala k ) + h k (R,ty(a k al)\, (10) 

k 



where h(R,t) = (2g/ui k ) e 
Tj — r m . If for all j 7^ m, 
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sin 2 (^Wfci) and R = 
: 0, then wc recover the 



U (t) = cxp(— iHst) cxp iH e ff(t) expG5s(t) , (4) independent noise condition Eq. (2). However, this is not 
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generally true. Although the phase factor e~ tk ' R tends 
to cause lots of cancellations, the remaining term tends 
to suppress this effect. Recall that sin 2 (w k t/2) jio\ — > 
(tt/2) 5 (u k ) t, for t > t m ax = Taax.{u)- 1 ,R/c,h/k B T} 
being the longest time scale. Roughly speaking, con- 
structive interference occurs (i.e., | (4>j4>m) I ~ I i&jrf'j) l> 
j ^ m) whenever the waves of disturbance has enough 
time to travel from spin j to spin to, i.e., t > £™ ax = 
max{|rj — r m |/c}. In the case of phonon bath, if the 
spins are located in a region of area 100/im 2 , then t™ ax w 
10 _7 s, which suggests that to maintain the independent 
noise approximation, the whole error correction proce- 
dure has to be executed at least 10 7 times per second. 

From the point of view of current experimental situa- 
tion, it is not unreasonable to consider the "worst case" 
scenario where we set 



4>m4>m) (11) 

for all j and to, and consider the impacts to quantum 
error correction. In this case, the probability (upper 
bound) P n < A 1 n of n-qubit errors would be enhanced 
(relative to the independent case P n = P") by a fac- 
tor {2n)\/2 n n\ w ^/2(2n/e) n for large n (e.g. n = 2 k 
for fe-level concatenated distance-3 codes). Consider the 
standard error threshold analysis based on the indepen- 
dent noise model. The failure probability is given by 
Pfaii = Pth(Pi/Pth) n , where P th is a threshold error 
rate, which depends on the coding methods and circuit 
designs. When Pi < P^, Pf a u can be made arbitrar- 
ily small by increasing the levels of concatenation. Now 
in the presence of the constructive interference effects, 
Pf -► \/2(2n/e)"Pf , hence 

■^•^(ts)' ■ (12) 

The concatenation method becomes inefficient when 
Pi ( e /^ n ) Pth and breaks down beyond that. Therefore, 
the interference effect imposes a more stringent threshold 
for quantum error correction. 

The good news is that, in the same limit, provided 
that the effective interaction [cf. Eq. (7)] is properly 
compensated, the correlated errors can be made vanished 
within the decoherence free subspacc (DFS) [12, 14]. In 
the ordinary analysis of DFS, one requires H$b |DFS) = 
U | DFS), where U is either zero or some unitary opera- 
tors acting only on the qubits. This condition, requiring 
zero qubit-qubit separation, is not assumed in our case. 
However, the qubits do effectively "see" the same envi- 
ronment. Since when the long wavelength modes k — > 
become dominant, the spatial qubit separations cannot 
be "resolved" . To justify this, we consider the subspace 
of states satisfying z j \DFS) — 0, e.g. all sym- 

metrized states with equal numbers of '0' and '1'. From 
Eq. (11), \P%b) dfs = ^ an< ^ hence by invoking the 
Bloch identity again 



It is now clear from Eq. 
dynamics is decoupled i.e. 



(1) that the spin-bath 
U{t)\DFS) <g> \Env) 



e iH eft {t) \ D FS) ® e~ lHBt \Env). Therefore, the idea of 
isolating the effective interaction gives us new insights 
about the method of DFS (answer to question (2)) and 
makes it physically more applicable. 

Case II: Bit-Flipping — In the following, we shall con- 
sider the effects of including Hs ^ 0. It is obvious that we 
reach the same conclusions if [Hs,Hsb] = 0. The sim- 
plest non-trivial case seems to be Hs = (A/2)Xj. 
As mentioned before, the coupling with the bath causes 
transitions between the two eigenstates of X, and is most 
efficient for Tiutk ~ A, as A sets the energy and hence 
the distance scale A* = Tic/ A for our problem. For 
weak spin-bath coupling Hsb, we employ the standard 
method of canonical transformation: H = e~ s He s = 
H S + H B + (1/2) [H S b,S] + where S satisfies the re- 
lation H SB + [H s + H B ,S} = 0. Define Lf = \±) (t\ = 



(Zj =F iYj) j2 and S = J2?=i Sj, we write (h = 1) 



Ofc 



(14) 



where T 3 {oj k ) = Lj / (A - u k ) - L+ / (A + u k ). The 
evolution operator can be expressed as U (t) = 
expS'exp(— iHt) exp (— S). The (lowest order) effec- 
tive interaction can be extracted from the commutator 
[Hsb, S], it turns out to be of the ZZ form: 



Ha. 



ZjZ„ 



k 



2g 2 



LJ k 



cos [k ■ R 



(15) 



which reduces to the first term of Eq. (7) when A — ► 0. 
The second term, containing sin(u;fct), emerges when we 
combine and rearrange the operators in U and transform 
it into the form (with Hb — ► Hb + Hs) in Eq. (4). Cru- 



cially, since all diagonal matrix elements (n\ Hj 







are zero in the eigenbasis of Hs, the term Hj m does not 
grow linearly in time (bounded by a factor 1/A) and be- 
comes negligible [15] for At/fr^> 1. 

We shall now consider the interference effects. For 
the moment, we neglect the higher order correc- 
tions and any effective interaction generated, and 
write U{t) w e ~*(Hs+H B )t e S(t) e -S(o)^ where s ^ = 
e i(H s +H B )t S (Qj e -i(H s +HB)t i We are mos tly interested 
in the limit where the energy non-conserving terms in 
Eq. (14) becomes relatively small. Then, the counter- 
part of Eq. (8) is 

N 

exp G SB (t) = I] CX P [ iY iV V j (*)] CX P \ Z M (*)] . ( 16 ) 

3=1 

where, with r\ = for and r\ = 1 for (p v -, 



(expG S B) DFS = 1 



(13) 



r 



■ V (t)^^[/ fc (r„t)4-(-l)"A(^,tr 



a k 



(17) 
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and f k ( rj ,t) = [g/{uj k - A)]e- ?;£ ^[1 - e^K-A)*]. Now 
the same argument, about the error upper bound, follow- 
ing Eq. (8) should go through the same way. However, 
the two-point correlation functions {(pj' y (p^}, which 
cause spatial error correlation, vanishes when the qubits 
are separated sufficiently far apart At s /% 3> 1, where 
t s = \fj —r m \/c. For example, (fj<Pm) is almost the 
same as in Eq.(lO), except the replacement hk (R,t) — > 
hk (R, t) where 

h k (R, t) = - 4g ' 2 e- itn sin 2 [\ (oj k - A) t] , (18) 
[ujk - A) 

When t i max = maxj^r 1 , t s , ft/fe^T}, we have 
(fjfm) ~ (PjtPj) sm (At s ) / At s . In fact, to have the in- 
terference terms get cancelled, the condition At s /h 3> 1 
could be too strong. It is likely that for At s /h~ ~ ir, the 
correlation functions start to change signs and cause lots 
of cancelations. This however should depend on the spa- 
tial distribution of the physical qubits. Under these con- 
ditions, we expect that the independent noise approxima- 
tion is valid (answer to question (1)) for the spin-boson 
decoherence, provided that the higher order terms, so far 
neglected, are also small compared with that generated 
by the single-error terms discussed [cf. Eq. (14)] above. 
Below we shall give a heuristic argument to justify that. 

In the canonical transformation, higher order correc- 
tions terms are generated from the series of commutators 
[[[-ffss; S] , S] , 5]. For each term, we denote (n, m) to 
represent the number n of Pauli matrices (or 'weight'), 
and m the weight of eft and a, for example (1, 1) for both 
S and Hsb, which are 0(g). Taking the commutator 
with S, each term changes from (n, m) to cither (n, m+1) 
or (n+£, m— 1) where £ = {0, ±1}, corresponding to tak- 
ing the commutator for the Pauli matrices and the bath 
operators respectively. The case (n+1, m— 1) is the only 
case causing an increase in the weight of the Pauli ma- 
trices, with a cost of increasing one power in g (making 
it smaller). The effective interaction (2,0) in Eq. (15) is 
an example of this. Starting from Hsb, i.e., (1,1), apart 
from (2,0), any (n > l,m) term is at least 0(g n+1 ), 
which is smaller than the contribution obtained by S n 
which is 0(g n ). We therefore conclude that the higher 
order corrections are negligible for sufficiently small g. 

Lastly, we remark that direct qubit-qubit interaction 
can also be included in H§ following the same line of 
thought, and the results should be qualitatively the same, 
except (not surprisingly) that the errors on the two in- 
teracting qubits become strongly correlated. On the 
other hand, it can also be generalized to determine the 
non-Markovian behaviors of the bath between the error- 
correcting cycles. In this case, the space-time correlation 
functions, (tpj (t a ) <p m (h)) where t a ^ are needed to 
be considered. Details to be given elsewhere. 

Conclusions — In conclusion, we have illustrated the 



general structure and essential features of the spin-boson 
decoherence for quantum error correction. The noise cor- 
relations are determined by the relevant two-point corre- 
lation functions. In particular, we have considered two 
special cases where the correlated errors can either be 
generated freely in the dephasing case or be suppressed 
completely (when At s /h 3> 1) in the bit-flipping case. In 
the former case the error threshold becomes much more 
stringent [cf. Eq. (12)] , but it could be handled with a 
modified decoherence free subspace (DFS) method. In re- 
ality, environment induced decoherence may exhibit fea- 
tures of both cases. It is likely that the combination [16] 
of the two methods, namely strong local field (hardware) 
and error-correcting codes together with DFS (software) , 
could help to minimize more general errors. Based on 
the results of this work, spatially correlated errors do 
not seem to be a fundamental barrier for fault-tolerant 
quantum computation. 
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